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Abstract 

The universal method of expansion of integrals is suggested. It allows in particular to 
derive the threshold expansion of Feynman integrals. 
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During the last two decades the efficient techniques were developed [Q]-^] for expansions 
of Feynman integrals in external parameters. These techniques allow to perform practical 
analytic calculations in quantum field theory when exact integrations are not possible. 

Recently the prescription was suggested for expansions of Feynman integrals near 
thresholds but its consistent derivation is missing, see also [Q. 

In the present paper we suggest the universal method of expansion of integrals - 'the 
method of direct expansions'. In particular it allows to derive consistently the threshold 
expansion. 



To regularize integrals we will use dimensional regularization |10| for convenience (but the 
method is regularization independent). The dimension of the momentum space is D = 4 — 2e 
where e is the parameter defining the deviation of the dimension from its physical value 4. 

To demonstrate the essence of the method of direct expansions we begin with the following 
simple integral depending on some external parameter t 
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where e is the parameter of dimensional regularization. 

We want to expand this integral in small t before the integration. Let us directly expand 
the integrand in small t 
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_i + t_(l + e)i- + 0(t3). 

In the absence of dimensional regularization, i.e. for e = 0, this expansion would produce 
non-integrable singularities at x = 0. But dimensional regularization makes these singular- 
ities integrable and this direct expansion correctly reproduces the small t expansion of the 
exact result in the right hand side of eq.(||). Only the term is invisible for the small t 
expansion (i.e. its expansion gives zero in any order) due to the basic property of dimensional 
regularization 

= for t = or t = 00, n / -e, (3) 

where n can be an arbitrary number. 

One can say that the term is a non-perturbative term for the small t expansion. Thus 
the expansion (|2|) is correct but not complete, since it has the non-perturbative term which 
can be not small for small t, e.g. its expansion in e gives ln{t). 

Let us now expand the integral Ii in large t. We again directly expand the integrand. 
This means that one factorizes the factor and then expands the integrand in a Taylor 
series in 1/t or, in other words, one expands the integrand in small x 
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Again this expansion correctly reproduces the expansion in large t of the exact result in the 
r.h.s. of @). The expansion in large t for the exact result means that one first factorizes the 
non-perturbative factor ^ and then expand the rest in the Taylor series in 1/t. The expansion 
@) is complete, i.e. there are no non-perturbative terms for it. This is because the terms of 
this expansion are integrable even for e = 0, i.e. in the absence of dimensional regularization. 

Let us generalize the above considerations. This simple example reflects the principal 
property of expansions of integrals: an expansion of a given integrand produces a correct 
result if all terms of this expansion are integrable. (For the integral Ii it is dimensional 
regularization which makes all terms of expansions integrable. For some integrals it can 
happen with the help of other regularizations, e.g. analytic regularization.) A correct result 
means that it correctly reproduces a corresponding part of an exact result for a given integral. 
But the essential point is that a correct result does not always mean a complete result. There 
can be non-perturbative terms in an exact result which are invisible for some expansions, i.e. 
expansions of these terms give zero in any order. (In the above example it was the term 
which was invisible in any order of the expansion in small t.) To get a complete expansion of 
a given integral one should construct a complete set of direct expansions of an integrand to 
reproduce all necessary terms. 

To see how the method works practically let us apply it to a Feynman integral. We will 
consider the following two loop propagator Feynman integral with an external momentum q 
and a mass M 

J {k- q)^k^[{k - ly - M2](g - /)2/2 ' 

where integrations over loop momenta k and / are performed in D-dimensions. Here and 
below in the paper we omit the "casual" iO in propagators for brevity. We will expand this 
integral in large M, i.e. in small q^ /M"^ . 

Let us first obtain the part of the complete expansion proportional to l/(g^)^'^. Then we 
will construct also the direct expansions for the remaining parts proportional to the factors 
l/{M'^q^Y and 1/(M^)^'^ allowed by power counting of non-integer dimensions due to non-zero 
e. One should construct direct expansions for all parts with different non-integer dimensions 
of external parameters. This is the essence of the method. 

To get the part proportional to l/(g^)^'^ we expand the integrand directly in large M. Or, 
in other words, we expand the massive propagator in the Taylor series in momenta k and I 

-l^ {k-lf^ 
M2^^ M2J ■ 



{k - If - M2 



Then we get 



(7) 



'ik-qyk^{q-iyP [{k - If - M-^ 

This direct expansion generates ultraviolet divergences, but dimensional regularization makes 
them integrable. Hence this is the correct expansion, i.e. it reproduces the corresponding 
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expansion of the exact result in large M (one first factorizes l/M"^ in the exact result and 
then expands the rest in small /M"^). By power counting of the non-integer dimensions due 
to non-zero e one finds that the integrals in (j^ are indeed proportional to l/(?^)^'^ since they 
depend only on one external parameter q. 

The complete expansion for the original integral I2 contains also terms invisible for the 
direct expansion (0). In particular the part proportional to 1/(M^)^ is invisible in this ex- 
pansion because of the basic property of dimensional regularization (^) which in this concrete 
case reads l/(iVf-^)'^+" = for M = 00. To get this part one can make the substitution 

= MK^ in the original integral I2 where is some new momentum. Then the invisible 
for the large M expansion factor l/{M'^y is explicitly factorized by the integration measure 
d^k = M^d^K and the integral takes the form 



d^K d^l 



(M2)3 J {K- q/MyK^[{K - l/Mf - l](g - IfP ' 
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In this expression one can directly expand the integrand in large M getting the necessary 
contribution proportional to 1/{M'^Y. (Then one can return to old variables making the 
substitution K^^ = k^/M.) One can see that the same direct expansion can be obtained 
from the original integral (^) by expanding the propagators carrying the momentum k in the 
Taylor series in momenta q and /. Thus one gets the second contribution 
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where Tq and T; are the operators generating the Taylor series in q and / correspondingly. 
Here the integrals over k and / factorize. The integrals over k depend only on one external 
parameter M and are proportional to 1/{M'^Y. The integrals over I are proportional to 
l/{q^Y. Thus we get the desired contribution proportional to l/^Ad'^q^y. One should also 
take into account that there is a symmetrical contribution of this type. It can be obtained 

from the original integral ^ if one expands the propagators carrying the momentum I in the 

(2) 

Taylor series in q and k. So the contribution I2 comes with the factor 2 to the complete 
expansion. 

(2) 

One more way to get the same contribution is via the substitutions q^ = and 
k'^ = A-fC^, where A is an arbitrary parameter. Hence the factor l/A^*^ is explicitely factorized 
by the integration measure d^k = X^d^K. Then one can directly expand the integrand 
in A and see that remaining integrals are proportional to the necessary factor 1/(AQ)2^. 
Substitutions of this type can be useful in expansions of more complicated integrals. 

We have obtained two contributions (0) and (|9|) to the complete expansion of I2. But 
the part proportional to 1/(M2)2'^ is still invisible in these equations. To get this part let us 
expand the integrand of the original integral (|^) directly in the Taylor series in the external 
momentum q. Then we get the third contribution 



(3) 



d^k d^l 



(10) 



4 



It is clear from power counting of non-integer dimensions due to non-zero e that integrals in 
the above equation are proportional to the desired factor 1/(M2)2'^ since they depend only 
on one external parameter M. 

Summing three contributions (0), (^) and (|To| ) we get the resulting expansion 
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One can ask if there are other non-perturbative terms in the complete expansion for the 
original integral I2 which are invisible in the above direct expansions, e.g. terms proportional 
to {q'^Y / {M"^)^^ . Such non-minimal terms (i.e. terms containing factors 1/(M2)^'^ or l/{q^Y'' 
where / is larger then the number of loops of the integral) could be present by power counting. 
But if one tries to construct a direct expansion of I2 which would generate such terms then the 
result is zero, which means that such terms are absent. In the case of the threshold expansion 
the terms of this non-minimal type appear for two loop propagator diagrams. 

The result (|l^) agrees with the prescription for the large-mass expansion of Feynman 
integrals explicitely formulated in |5|,[^. The essential point of the present paper is the 
compact derivation of this expansion. 

We will apply now the method to the case of the threshold expansion of Feynman integrals. 
Let us consider a scalar one loop vertex integral containing two massive propagators with a 
mass m and one massless propagator. This integral appears e.g. in the form factor for the 
decay of the virtual photon into a quark-antiquark pair ^*{q) Q{pi)Q{P2) with momenta 
vi = = ''^^j q = Pi + P2- One defines also the relative momentum p = {pi — p2)/2 with 
p ■ q = 0. The threshold region is characterized by 
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and one considers the expansion in small y, i.e. in small y/q'^- The corresponding integral is 

d^k 



[{k + 9/2)2 _ ^2]p _ g/2)2 _ ^2](^ _ py 
d^k 



(13) 



(A;2 -\- qk — y){k'^ — qk — y){k — p)2 

The expansion for this integral was obtained in Q but its consistent derivation was missing, 
see also §]. 

Below we give the derivation of this expansion using the method. One should again 
construct direct expansions for all possible terms with different non-integer dimensions of 
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external parameters. Let us first get the part proportional to l/{(f'Y. We directly expand 
the integrand in the Taylor series in y and p. All terms are integrable due to dimensional 
regularization and one correctly reproduces the expansion of the exact result for the integral 
I3 in small y. The leading term of this expansion is 

3 J {k^ + qk){k^ - qk)k^ W) 2e 1 + 2e " ^ ^ 

One can also calculate the higher order terms in y using e.g. Feynman parameters. 

But one can see by power counting that the complete result can contain also the part 
proportional to l/y*^. It is invisible for the small y expansion due to the basic property of 
dimensional regularization (^) which now reads l/y^"*"" = for y = 0. To get this part we 
should construct the direct expansion in large q. For this purpose it is convenient to go to 
the frame q = {qo, 0),p = {0,p)- Then the integral is 

' " I {k^ + q,k,-y){k^-qoko-y){k-py ^^^^ 
Let us make the substitution qoko = tq. The integral takes the form 

1 f drod^-^k , , 

h = — — z -2 z • (16) 

qoJ fy-k-^ + ro- y)(^ -k^-ro- y)C^ - (k - p}^) 
% % % 

We can now expand the integrand directly in large q^, i.e. in small Tq/^q. The resulting inte- 
grals depend only on one external parameter y producing thus the necessary part proportional 
to 1 /y*^ . The leading term of this expansion reads 

if) = 1 / '"^f"-'^^ . = (17) 

J {-k"^ + ro - y){-k'^ - To - y){-{k - p)"^) 

1 %/^r(e+l/2) 

Note that higher order terms in the expansion in large go give zero since positive powers 
of Tq produce (after integrations over rg) massless tadpoles which are zero in dimensional 
regularization. 

The exact result for the integral /s can be obtained by using e.g. Feynman parameters. 
Then one can see that the sum of the direct expansions (14) and ([T7| ) indeed reproduces the 
expansion of the exact result 

>h = -TT:^ 2Fi\ I + e, -;-—]= (18) 



i7r^/2 ' 2 \^2' ' '2' 

4y+^ 1 ^r(l + e + j) 1 / 4yV' 1 1 ^r(e + l/2 



f) 2e^^ l + 2e + 2j j! V y'+^Z^ 2e 
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In this way one obtains the threshold expansion by the method of direct expansions. 

The author is grateful to the collaborators of the Institute for theoretical particle physics 
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